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C\J ■ Abstract. The influence of higher order (stringly inspired) curvature corrections to the classical 

, General Relativity spherically symmetric solution is studied. In string gravity these curvature 

I corrections have a special form and can provide a singular contribution to the field equations 

' because they generate higher derivatives of metric functions multiplied by a small parameter. 

■ Analytically and numerically it is shown that sometimes in 4D string gravity the Schwarzschild 
CO , solution is not recovered when the string coupling constant vanishes and limited number of higher 

order curvature corrections is considered. 
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■ 1. Introduction 



The idea of higher order corrections to the Lagrangian of considered system became a 
rather common approach in modern theoretical physics. Sometimes, these additional 
corrections with higher order derivatives can drastically modify (through a singular 
contribution) the solutions of the corresponding Euler-Lagrange equations. This singular 
contribution does not disappear, even when higher order corrections vanish. This fact 
can be illustrated with a simple example. Let us consider the following Lagrangian: 

L = L{q)+e'l{q), (1) 

where 

m = f . 
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Let the initial conditions be 



g(0) = 

m = 
m = 



V, 







6q. 



As a result, one obtains two types of solutions. When /c = 1 it has the form 

q = Qo + Vot — Sq sm(t/e) 

and tends to the nonperturbed solution q = qo + Vot when £ — 0. When A; = 2 it looks 
like 

q = qo + Vot + 6q sh.{t/e) 

and has a singular contribution even when e ^ 0. 

One meets theories with higher order curvature corrections rather often. One of the 
first examples is the wellknown Wheeler- Feynman electrodynamics [|I|. It is based on 
the action that is parametrized in the form of two dimensional integral. Such integral 
can be approximated by the set of one dimensional integrals, so, an infinite set of time 
derivatives appears 0. Modified theories become unlocal and lead to the equations 
with the higher order derivatives up to the infinite orders. It is important that there 
is no problem of singular contribution from higher derivatives because the solutions 
are limited and are infinitely differentiable. Therefore, the lagrangian of infinite order 
approximates the lagrangian of second order which provides the regular contribution of 
higher order corrections. 

We would like to point out that the problems with new degrees of freedom, and, 
so, appearing of nonphysical "runaway" solutions are really actual for one special class 
of the models. There are the models where the corrections contain high but finite 
order of derivatives ||^, where the singular contribution of higher order correction is not 
compensated by the contribution from next order (s). The most natural approach to 
such class of models is to introduce additional limits on solutions availability and to 
consider only such ones that have a regular part (expanding in Taylor series) relatively 
the numerical parameter before the higher derivatives [Q. If one considers such condition 
as necessary one there would be no problem of higher order corrections (In Ref. [^] the 
problem is treated in such a manner). 

Here it is necessary to emphasis that higher order curvature corrections do not 
automatically cause singular unhmited solutions (see the example higher, Eq. (p). So, 
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it is desirable to study the conditions of appearance of new "runaway solutions" and their 
status in the four dimensional (4d) gravity models with higher curvature corrections 
caused by string theory more carefully. It is important to extract the properties of 
the solutions (in this paper we restrict our consideration with spherically symmetric 
ones) when only limited (finite) number of corrections are considered. Further, during 
consideration of real physical problems it would be better to restrict consideration with 
models where higher order corrections (that are necessary to make a little step from 
General Relativity boundaries of applicability) do not cause singular contribution to the 
regular part of the solution. The stability of the solution under the initial conditions on 
the infinity is strongly required from the physical grounds. 

In the perturbational approach, the 4d effective (string inspired) gravity action has 
the following form (we work in Planckian system of units where h = c = rripi = 1) 



where R is the Ricci scalar, is a dilatonic field, A is a string coupling constant (an 
expansion parameter which is proportional to string tension a' with the numerical 
coefficient depending upon the type of string gravity) and k {i = 2, 3, . . .) are higher 
order curvature corrections. In this paper we study the class of the models where higher 
order curvature corrections consist from Riemannian tensor R^uap products. Such type 
of models are widely discussed now from different aspects ||^. This study focuses on 
asymptotically flat spherically symmetric static black hole solutions. 

It should also be pointed out that the zeroth measure of the initial conditions 
manifold (only asymptotically flat, for example) is not very realistic, especially when 
new observational results in cosmology are taken into account. The highly probable 
existence of a global nonvanishing cosmological constant A requires the extension for 
the manifold of initial conditions. One can neglect the A influence on the global structure 
of the solution because it has zeroth order contribution relatively higher derivatives, so, 
it's influence is important only in the regions with rather small curvature. In black hole 
topologies where its influence is neglectible, it is possible to restrict the consideration 
by the requirement of stability (of Lyapunov type) under the initial conditions (and, 
hence, to take A existence into account) at the inflnity. 

The most convenient choice of four dimensional metric is, therefore. 



The solution of corresponding nonperturbed fleld equations is the wellknown 
Schwarzschild one 



S 




(2) 




(3) 



A 




r 



a = a\r) = l, (4) 
= 0°(r) = 00 = const, 

where M is black hole mass. 

When the correction I2 = RijkiR'^'^^ — ARijE}^ + in (j^) is taken into account 
(which is the wellknown 4d curvature invariant named as Gauss-Bonnet term), the 
corresponding solutions do not contain higher derivatives. It is necessary to emphasize 
that the discussed black hole solution behavior strongly differs from usual Schwarzschild 
one under the event horizon as it was shown in detail in Ref. 0. When the metric (|^) 
is substituted in (Q), the action becomes [§ 



S = dr 



Lrir, A, a, 0, A', 0', A) + X%ir, A, a, 0, A', a', A", A) 



(5) 



where T = b,h,s corresponds to the three string theory types considered in this paper: 
bosonic (b), heterotic (h) and superstring II (s). 

The basic Lagrangian does not contain higher derivatives. 



Lb = Lft = Lo(r,A,a,0,A',0') + ALGB(r,A,a,0,A',0'), 

AV + A + - Ar2(0')2 
Ls — Lq — , (bj 



a3 



Lq is simply the contribution resulting from General Relativity. 

The next order curvature corrections to /x have the following form |0: 

3 1 

+ R^\pR''%pR'^%,) + 0{X), 
k = C2v^e-6^A - liR^.o^pR^"''^? - -R^.^'R,s'"^Rpr,''^Ro.r 
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+ ■^^RpLV^^^RaP^R^ra^Rp^'^^ + RpLu'^'^RaP^^RprP^R'lS^'^) + 0{X), (7) 

Is = C3v^e-6'^A + 0(A), 
where 

^ = C{^)[Rpupr,R'^''^'^ {R'^" SpRa-i^"^ — '^R'^'^SaRfS^^'^)], 

Rfiu — RpyR^" 1 C(3) is Riemanian zeta-function, Ci^2,3 are numerical coefficients. 

The influence of higher order curvature corrections to the behavior of spherically 
symmetric static solutions was considered in where a perturbed solution close the 
nonperturbed one was studied. It was obtained by a specially developed and coded 
numerical iteration method. The only indefinite point of the perturbed solution occurred 
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near the particular point (internal black hole singularity of a new type at finite 
distance from the origin, provided by Gauss-Bonnet term, see Ref. 0). If one works 
in the frames of such an approach, the singular contribution disappears. From a pure 
mathematical point of view, the solution tends to some particular branch that begins 
from a particular manifold of initial conditions [IU| with null measure. Even small 
changes of initial conditions cause an appearance of a singular contribution that does 
not vanish when £ — 0. 

Finally, the main aim of this paper is to understand under which circumstances a 
singular contribution is induced in three different closed 4d low energy string gravity 
models on the Schwarzschild metric background when radial coordinate r is rather large, 
when the only finite number of higher order curvature corrections are taken into account 
(as the complete expansion is unknown). The paper is organized as follows: In Section 2 
we show some general mathematical considerations of higher order singular corrections 
behavior, in Section 3 we apply these results to bosonic, heterotic and superstring II 
(SUSY II) closed 4d low energy models. Section 4 contains the SUSY II numerical 
investigations. Section 5 is discussions and conclusions one. 



2. General theory of singular contributions 

Models with corrections including higher order derivatives can be generically written as: 

L{t, X, X, x) = L{t, X, x) + e'^l{t, x, x, x), (8) 

where e is an expansion parameter, t is real variable, x{t) is a smooth vector of n- 
dimensional manifold and L, I are smooth functions. 
The corresponding Euler-Lagrange equations are 

dL _ d^dL ^ d^dL _ ^ 
dx dt dx dt'^ dot 
One can rewrite these equations as follows 



X- 



dxdx dxdx 



+ $(t,x,i;) + 



-X 



(3). 



dH 



dxdxdx 



T7X^^^ + x^^^H(t, X, x, x) + \E'(t, X, X, x) 



(9) 



where and are symmetric matrixes of second order derivatives of L and I 

dxdx dxdx ^ 

(respectively) whereas $, S and \& are vector matrixes with second and third derivatives 
of L and I (respectively). 



Introducing the following change of variables: x = y, x 



z, ex 



(3) 



V and 
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considering the matrix as a nondegenerated one, we can write @ in the form 

j ev = F{t, X, y, z, e), ez = v, 
\ y = z, X = y. 

where 

F{t,x,y,z,e) = 



(10) 



dH 
dzdz 

dzdz 



n -1 



{t,x,y,z) 



z^:^^{t,x,y) + 



dydy 



dH 
dzdz 



-1 



dzdz 



dH ^ ^ ^ 
dzdzdz 



^i{t,x,y,z). 



dH 
dzdz 



-1 



v'~{t,x,y,z) 



When e = the system (^) has a degenerated solution X, 

X = y = y'^it), z = z^it), v = t;°(t) = 0, 

which corresponds to the nonperturbed one of Euler-Lagrange equations @ when the 
initial conditions are x^(to) = Xq, y^ito) = yQ. This means that the equations (p!0|) 



(11) 



belong to A.N.Tikhonov's class |]TT|. According to the related theorem, the closeness of 
perturbed and nonperturbed solutions can be obtained, with initial conditions in the 
neighborhood of xq, yo, x"{to), x"'{to), if all the real parts of eigenvalues u of matrix 



dz 



\I J 



are negative. This matrix is calculated on nonperturbed solutions X (|11]), I is identity 
n X n matrix. 

After some algebra, the equation on eigenvalues u can be represented as 

dH d^L 



det 



dzdz dydy 



0. 



Therefore, to provide the convergence of the perturbed solution to the nonperturbed 
one, it is possible to require only Re{uj) = 0. This corresponds to the case when the 
value 



u{t) < 0, 

where u{t) can be calculated from the following equation 
dH d^L 



det 



u{t) 



dxdx dxdx 



0. 



(12) 



(13) 



x=x{t) 
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We shall call u{t) the singular index. Eq. (p!3| ) is valid on the interval t G [to?^] where 
the solution of nonperturbed problem exists when e = 0. 

It is important to emphasize that the condition (|12|) is not enough to ensure the 
absence of any singular contribution to Eqs. (P), but when the condition (|I^) is broken 

a singular contribution of order exp^- Re^Ju{t)^ appears necessarily. 

Switching to the string gravity perturbed Lagrangian 

L = L{r, A, a, 0, A', 0') + l{r, A, a, 0, A', a, A") 

and taking into account that 0(r) is not perturbed, one can conclude that the condition 
(O) on u(r) can be fulfilled if 



\uP-QU- 



where 



and 



:AO(r), cr=aO{r), (/)=0O(r) 



P 



0, 

Ia"A" la'A" 



Q 



(14) 



3. Singular indexes of low energy effective string gravity 



To investigate the problem of singular contributions in string theory Lagrangians with 
higher order curvature corrections in spherically symmetric space times ,it is necessary, 
as a first step, to apply the formulas (|^) and to determine the accurate form of the 
main Lagrangian and higher order curvature corrections. The software packages from 
MAPLE and REDUCE were specially developed and used to obtain the corresponding 
singular indexes. With the help of Eq. (|14D and after the substitution of Schwarzschild 
values of metric functions in the limit r/M — oo, we get the following asymptotic 
formulas. 

For bosonic case 



u 



1,2 



For heterotic case 



u 



u 



576 



0.008773e^'^°r^ 



-0.022916^"^^^ 



r / r 

ho — 

M \M 



J.2 ( \ 

M2 ^""Im^ 
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For SUSY II case 




Considering these asymptotical formulae, one can conclude that only SUSY II 
does not induce any singular contribution to the solution in the neighborhood of the 
Schwarzschild one , even for big r. To understand the perturbed solution behavior 
in SUSY II case, it is necessary to perform a numerical investigation of this solution 
dependence upon different sets of initial conditions A"(ro), A"'(ro), cr'(ro), (T"(ro), that 
are close to the corresponding Schwarzschild values AQ(r) = — 4M/r^, AQ"(r) = 12M/r^, 



4. Results of numerical investigation of perturbed solutions in SUSY II 

To investigate the difference between the perturbed SUSY II solution and the 
nonperturbed Schwarzschild one , we solved Eqs. as the Cauchy problem starting 
from initial point tq using initial conditions from accurate Schwarzschild solution. A 
7th order Runge-Kutta code was especially chosen for the integration of the resulting 
system. Some computing problems occurred in the range of big tq/M. For instance, 
when A = 0.1, M = 1 for tq > 6 (we work in Planckian units) the system becomes so 
rigid that it can not be solved. We could integrate these equations only in the range 
that was situated up to the event horizon at r/M e [2.1, 5.2] and for M < 100, because 
for higher mass values and higher r/M values local influence of the corrections became 
so small that calculations stopped at the first integration step. 

For different values of rg, the difference of A(r), ij{r) and 0(r) from corresponding 
Schwarzschild values with high accuracy (not less than integration method one) were 
proportional to (r — Tq)^, so. 



where C is numerical coefficient. 

We investigated the dependence of C as a function of the initial integration point tq 
and mass M when A is fixed and equal to O.lf. So, the behavior of C(ro, M) is described 
by the following fit 



<{r) = <(r) = 0. 



Sx = C{r — TqY, 



(15) 




f The increasing of M when A is fixed is equivalent (up to dimensionality) to decreasing of A when M 
is fixed 
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It is important to emphasize that I3{M) > 10 has very weak dependence upon M. For 
instance, the results for 0(r) {(po{r) 7^ 0) are presented in the Table 1. 



M 


a 




1 


134. 


13.53 


2 


0.411 


13.59 


15 


4.39 ■ 10-8 


13.81 


50 


3.29 ■ 10-12 


13.81 


100 


1.49 ■ 10-1^ 


13.92 



Table 1. The dependence of the numerical fit coefficients of expansion a and /? versus mass M 



If one extrapolates formula (|T^) for 6x (difference between perturbed SUSY II 
solution and nonperturbed Schwarzschild one) in the complete range of r/M G (2, cxd), 
after obtaining the asymptotically flat solution (starting from tq 00), one would 
conclude (for the fixed r and (3 from Table 1) that: 



lim C(ro, M)(r - roY = lim C(ro, M)rf. 1 + o — = 0. 



Then it is possible to estimate 5x{r) as 



2 4«(M) f (3 r\-\ 



For example, at the event horizon r = rh = 2M : 

6x{2M) < D{(3)a{M)M\ 

where 



(13 - 2)2 \I3 - 2 J ■ 

If the value of parameter (3 does not differ apparently from those presented in Table 1 
and the dependence a(M) is approximately the same, one can conclude that 

lim a{M)M^ = 0. 

Af— ♦00 

So, the value of the perturbation becomes zero when relative weight of higher order 
corrections vanishes. 

Finally, this, as it seems to the authors of the paper, means that the 4th order 
curvature correction in SUSY II does not provide any appreciable differences from the 
Schwarzschild solution, at least up to the event horizon. 
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5. Discussion and conclusions 

Based on our preliminary results in the framework of models with higher order curvature 
corrections made of pure Riemanian tensor products, it is reasonable to assume that 
the only model that does not provide any singular contribution is the SUSY II one. As 
curvature corrections of next orders do not produce new higher derivatives, in a case if 
singular contributions are absent in the 4th perturbation order, they would be absent 
in all perturbative orders. It is impossible to draw more accurate conclusions in this 
framework. Our results show the existence of singular contributions in bosonic and 
heterotic realizations of string gravity and give hope that SUSY II model can be free 
from singular contributions, at least for spherically symmetric static space times. 

It is necessary to emphasis that the loop expansion of string gravity (the origin of 
higher order curvature corrections) initially consists from the terms with increasing order 
of derivatives |T^. But working in the frames only of one considered approximation 



(for example, 3rd, 4th, ...) it is possible to make some manipulations taking usage 

from the fact that corrections contain the products with dilaton. So, analogously to the 
transformation that was used when the Gauss-Bonnet term by itself was changed to its 
current [j^, one can make partial integration of the correction according to the scheme: 

I e"*Vni?" = (e"*V„_ii?") - I nV$e"*Vn-ii?". (16) 

Making this procedure "n" times one can shift differentiation from Riemannian tensor to 
the dilaton and finally result with the terms like R^j^i without derivatives. Then, as one 
is working in the limits of corresponding correction, the term with dilaton derivatives 
can be transferred to a next order as it is shown in Formula (10) at Ref. |TB[. Hence, 
in the frames of chosen order the expression becomes free from derivatives. 

So, the result of our study is that when the limited number of higher order curvature 
corrections are taken into account a spherically symmetric solution can disappear. 
Finally, the singular contribution of the higher order curvature corrections in the low 
energy limit of string gravity, in the frames of Lagrange approach, represents a problem. 
Anyway, as the consideration of higher order corrections is necessary to extend the 
boundaries of applicability of classical gravity (and to obtain preliminary directions on 
some effects, that have quantum gravity nature) this can be avoided by two possible 
ways: 

• to suggest additional principles for choosing the special initial conditions for 
the higher derivatives, providing the closeness of the perturbed solution to the 
nonperturbed one (at least for big r) that include all the physically interesting 
cases (asymptotically fiat, dS/AdS); 

• not to restrict the consideration by the frames of chosen perturbative order and 
consider the "rest term" of series expansion (this means that the application of 
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Formula (10) from [|I3] and all the same ones are not directly allowed in black hole 
cases) . 
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